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Abstract: A new class of neutrino dark energy models is presented. The new models are 
characterized by the lack of exotic particles or couplings that violate the standard model 
symmetry. It is shown that these models lead to several concrete predictions for the dark 
energy equation of state, as well as possible effects on the cosmic structure formation. 
These predictions, can be verified (or disproved) with future experiments. At this point, 
the strongest constraints on these models are obtained from big bang nucleosynthesis, and 
lead to new bounds on the mass of the lightest neutrino. 
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1. Introduction 

It has been over a decade since observations of supernovae revealed that the expansion of 
the universe seems to be accelerating [1]. Today in addition to supernova observations, 
the acceleration of the universe, is supported by data from WMAP [2], SDSS [3] and more 
[4]. The source of this acceleration, the so called "dark energy" (DE), remains a mystery, 
since it cannot be any of the well known cosmological components. The only previously 
known energy source that can play the role of DE is the vacuum (or zero-point) energy. 
The magnitude of the zero-point energy density depends on the QFT cutoff. However, 
any natural cutoff, generates an enormous vacuum energy, up to 120 orders of magnitude 
greater than the value needed to account for the observed acceleration of the universe. 
Thus, in addition to the strange nature of DE, one must account for its "tiny" value 
(Pde ~ (l0 -3 ey) 4 ), the same value responsible for the "cosmic coincidence" - the fact 
that today, DE is of the same order of magnitude as the cold dark matter density, despite 
the fact that the ratio between them changes as the cosmic scale factor cubed. 

Countless theories have been presented to explain the source of cosmic acceleration and 
its value. Most of them include either modified gravity [5, 6] or exotic sources of energy 
[7, 8], deviating far from the standard model (SM). Since strong deviations from known 
physics are usually accompanied by numerous degrees of freedom, testing for these theories 
is very difficult, in spite of the advancing observational constraints on DE [2, 3, 4, 9]. 
In addition, most theories seem to be unable to naturally generate the correct value for 
DE, usually leaving it as another degree of freedom, resulting in the need for fine-tuning. 
Finding ourselves in such a vast landscape of alternative theories, emphasizes the need for 
a theory with minimal deviations from known physics and preferably no fine tuning. 
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neutrino masses are the closest known fundamental energy scale to the scale of DE 
[10, 11]. This implies that neutrinos might be connected to DE, and that a theory with an 
underlying link between them would have less additional degrees of freedom. 

This idea is not new, and various methods to link neutrino physics to DE, have been 
investigated [12, 13, 14, 15, 16, 17]. One of the leading models for neutrino DE is mass 
varying neutrinos (MAVANs) [18, 19]. This class of models includes coupling between 
neutrinos and a scalar field (the Acceleron). The DE "tracks" the matter density to 
a critical point, when the temperature of the cosmic background and the mass of the 
neutrino are comparable, and then becomes constant and generates the acceleration we 
witness today. All the models, mentioned above, have specific issues, but one issue is 
common to all of them: It is unnatural for the neutrino to be the only source of DE. For 
example, in the case of MAVANs, there is no reason for the Acceleron to couple only to the 
neutrino. Being a neutral scalar particle, coupling to a neutrino must imply coupling to 
the electron (and by the same token, the muon and the tau). A natural model of neutrino 
DE cannot include such an exotic coupling that distinguishes the neutrino from the other 
leptons, thus, breaking the SM symmetry. 

2. Constructing the Model 

In order to construct a natural neutrino DE (NNDE) model, it is essential to find a quality 
that separates the neutrino from the other leptons. The obvious would be the lack of 
electric charge. Models, in which the DE component is coupled to Majoranna mass terms, 
naturally account for the neutrino being the only source of DE. However, since Majoranna 
masses are typically very high, we lose the correlation between the mass and the DE scale 
and thus, the entire motivation behind neutrino DE. Therefore, it is necessary to find a 
different quality that would distinguish the neutrino. 

The three neutrino species, appear to be the lightest fermions. Particularly, the lightest 
neutrino is the only known massive particle that may still be relativistic in the cosmic 
background. The heavier neutrinos are not relativistic (based on neutrino oscillations 
[10]), but they are much closer than any other known fermion. Thus, we can base a NNDE 
model on the idea that relativistic particles contribute to DE, while non-relativistic ones 
do not. We construct a model, such that the energy density is 

Pde(«)= J2 U & (°)>Xi (°)>-) • M 

i=particles 

Here £i(a) , Xi( a ) ? ••• are scalar functions of certain parameters of the cosmic ensemble (re- 
ferring to the particle species i). £i(a) and Xi( a ) are almost constant in the ultra-relativistic 
regime and diverge as the kinetic energy of the particles drops to zero (alternatively, as the 
cosmic scale factor a —> oc). The function U must decay to zero as £(a),x(a), ... diverge, 
so there is no contribution to DE from non-relativistic particles. Also, to correlate DE to 
particle masses, U ~ m 4 in the ultra-relativistic regime, where m is the particle mass. In 
eq. (2.1), it was not specified, whether we sum over all elementary particles, fermions or 
only leptons. This is to be determined in a specific model and is irrelevant at this point. 
Eq. (2.1) is not restricted to a specific model and guaranties the following: 
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• No exotic couplings that violate particle physics symmetries are necessary, since DE 
couples evenly to every particle. 

• At any given time, DE is dominated by the heaviest relativistic particles. Heavier, 
non relativistic particles, have no contribution. Therefore, 

— Today, DE is dominated by the (lightest?) neutrino. 

— In the past, heavier particles contributed as well and the DE density was greater. 
This, will be referred to as primordial DE. 

The form of eq. (2.1) is therefore sufficient to explain why the neutrino is the only one to 
contribute to DE today. 

To proceed, it is necessary to find all the possible parameters £(a), x(a),... that satisfy 
the above mentioned conditions. The mean energy E must be a factor in all the parame- 
ters, since it defines the relativistic and non-relativistic regimes. E scales with expansion 
oppositely to what we have defined; it is constant in the non relativistic regime and diverges 
in the ultra-relativistic one. Therefore the energy must be suppressed by a quantity that 
constantly grows with expansion. Thus, two simple parameters come to mind, £ = El and 
X = E/T. Where / is the mean distance between the particles and T is the temperature 
of the cosmic background (T is not necessarily the CMB temperature. For example, for 
neutrinos T = T CM B (4/H) 1/3 ) 



3. NNDE Phenomenology 

At this point, we can find an expression for the DE equation of state (EOS), 

1 + wde o — • (3-1) 

3 Pde ocl 

Assuming, only a single neutrino species gives a significant contribution, substituting 
eq. (2.1) into eq. (3.1), we have 

i a fouae dud x \ 

For cosmic neutrinos [11, 17], we can write 

47Tp 2 dp \J m 2 + p 2 



/4:7rp dp 



( 27r ) 3 l + exp(p/r) 
J (2tt) 3 l + exp(p/T) 

(3.4) 



f 4irp 2 dp 1 \ 

J (2tt) 3 l + exp(p/r)J 



1/3 • 
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Model 


w a (w ) 


NNDE 


w & -2(1 + wo) 


basic MaVaNs 


w a = -3(1 + w ) 2 


Vacuum fluctuations 


= -3(1 + ^ ) 2 


Chaplygin gas 


w a = 3(l + w ) (l-(l + w ) 2 ) 


DGP 


w a = 3(1 + w ) (l - 2(1 + wo) (l + y/-w (l + w )) ) 


Quintessence Thawing 


-3(l + ^ )<^a<-(l + ^0) 


Quintessence Freezing 


-0.2w (1 + w ) < w a < -3^ (1 + w ) 



Table 1: The linear evolution approximation of the DE EOS, for different models. From top to 
bottom: NNDE, a simple MaVaN model [18], Vaccum energy model [15], Chaplygin gas [8], DGP 
[5], quintessence thawing and quintessence freezing [7]. 



Expanding the square root in eq. (3.3) via \J m 2 + p 2 — p J rm 2 /2p and substituting it 
along with eq. (3.4) into eq. (3.2), we have 

^-^^ir^^jt) ' (3 - 5) 



where, 

* - h {ymT * °- 268 and h = WW) * a152 

are numerical coefficients. U and its derivatives are functions of £ and %, both of which are 
constant to first order, in the ultra-relativistic regime. Therefore, in the ultra-relativistic 
regime T ~ a -1 is the fastest varying quantity in eq. (3.5), it follows that to first order, 
1 + ^de ~ o?. Expanding eq. (3.5) via powers of 1 — a we have, 

^de = m + w a (l - a) + w 2 (l - a) 2 , (3.6) 

with 

w a = -2w 2 = -2(l + w ) , (3.7) 

and 

, / dlnU n dluU\ fm\ 2 

Where To is the temperature at the present time. Eq. (3.6) is valid even if u>de deviates 
strongly from -1, as long as the ratio m/T is small, and thus the particles are relativistic. 

Verifying the relation w a (wo) is a good, parameter independent, method to test a 
variety of DE models. Table 1 shows the linear evolution approximation of the DE EOS, 
for different models. Some of these models are plotted in figure 1. One can clearly see the 
distinctive behavior of the different models. 

As shown, NNDE corresponds to the thawing range of parameters for w > — 1. A way 
to distinguish the two may be to analyze observational data, using eq. (3.6) rather than its 
linear approximation. 
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Figure 1: The linear evolution approximation of the DE EOS, for different models. The plotted 
contours signify the observational constraints, taken from [2]. Within the observational constraints, 
the different models are clearly distinguishable, aside from the DGP and Chaplygin gas curves. 

Another important prediction of NNDE models lies within its initial formulation. 
There is no natural way for neutrino to induce DE, without heavier leptons inducing it 
as well. We have eliminated the contribution of heavier leptons by constructing a model 
where only relativistic particles have significant contributions. However, this means that 
in the early universe, heavier leptons did contribute (primordial DE). The further we go 
back in time, the more contributions we have. In order for us to discuss the cosmic history 
in the presence of NNDE, we must choose a specific model, that is to determine a specific 
function J7(£, x). As a toy model, let us use 

17 = Cm 4 exp(-£), (3-9) 

where, C is a constant. I emphasize that the toy model is used merely for exact calculations 
and all the conclusions below are model independent. 

Since, unlike neutrinos, many fermions stay coupled to the cosmic plasma in the era 
where T ~ m (which is the era where the DE of the specific particle dominates the ex- 
pansion), the exact contribution of those fermions is slightly different. Figures 2 and 3 
show exp(— £) and exp(— x) as a function of m/T for coupled and decoupled fermions. It 
is visible that the general behavior of the model for coupled fermions, remains the same. 

Figure 4 shows the expansion history, in terms of the Hubble plot, for ACDM vs 
NNDE. On the Hubble plot, the mass hierarchy of the different leptons is clearly visible. 
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Figure 2: exp(— £) vs m/T for couples and Figure 3: exp(— x) vs m/T for couples and 
decoupled fermions. decoupled fermions. 



The "staircase" form of the Hubble plot is in correlation with the masses of the different 
leptons 1 . The first DE "step" (the one closest to the present epoch) is due to the fi neutrino, 
the next are due to the r neutrino, then the electron and so on. 



— NNDE 
— ACDM 




Figure 4: The Hubble plot for NNDE and ACDM models. 



lr rhe neutrino masses are unknown so this is a schematic plot. 



4. Constraining NNDE Using Big Bang Nucleosynthesis 



4.1 NNDE Constraints 

Other than using the relation dictated by Eq. (3.7), It is difficult (although not impossible) 
to constrain the NNDE model via usual astrophysical observations, including CMB. The 
reason is, that NNDE phenomenology does not differ much from other DE models, including 
the cosmological constant. The major difference lies in the primordial DE. However, since 
we do not know the neutrino masses, it is difficult to state, at which phase in the 'recent' 
expansion of the universe, do primordial DE effects show themselves. 

In fact, a well posed constraint on the NNDE comes from big bang nucleosynthesis 
(BBN). Since the BBN freeze-out occurs at a temperature of about ~ lMeV [20], it should 
be strongly affected by the electron-induced DE. Unlike the neutrino masses, we cannot 
maneuver our way around this issue. Both the BBN ratios [21] and electron mass are 
measured with good precision and a significant variation of either is out of the question. 

There are two ways to ensure that the presence of the electron would not have signifi- 
cant effects on the BBN, without abandoning the principles of NNDE: 

• The first way is to construct a model, where the contribution of any particle to DE 
decays before it becomes non relativistic (To be referred as type I NNDE). This would 
assure that by the time of the BBN freeze-out, the electron-induced DE is negligible. 

• Finally, it is possible to construct a model such that the contribution of any particle to 
DE, is suppressed by a factor of several orders of magnitude, so it is still proportional 
to m 4 , but lower (To be referred as type II NNDE). This can render the electron- 
induced DE low enough so it would not significantly affect the BBN. 

To proceed, we modify the toy model from Eq. (3.9) to 



Here, g is added as a parameter which tunes the energy scale, beyond which a particle no 
longer contributes to DE. 

4.2 Type I NNDE 

In this model, to assure that the electron-induced DE does not affect the BBN, we require g 
to be large enough, so that ex.p(—gEl) quickly decays for temperatures T < 3MeV. Using 
Eq. (3.5) obtained for the EOS, the above requirement can be written as 



At the same time, we require 1 + wq <C 1. Thus, the neutrino mass must be small enough, 
so that 



U — Cm 4 exp(— g£). 



(4.1) 




(4.2) 




3 • 10~ 5 eV, 



(4.3) 
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where T v — Tcmb(4/11) 1 ^ 3 . Figure 5, presents the constraints on the neutrino mass and 
the parameter g, obtained from BBN freeze-out (neutron to proton ratio) and the DE EOS. 
The BBN freeze-out ratio, was calculated numerically using a generalization of the model 
presented in [20], 

dy 

2y 2 (l + cosh(l/^)) X 
x exp (-A jT dy'H^y' 2 (y f + 0.25) 2 (l + exp(-l/y / ))) • (4.4) 

Here, H is the Hubble constant, A ~ 3.26sec _1 and y = T/Q, where, Q c± 1.29MeV is the 
difference between neutron and proton mass. This model gives a good approximation of 
the freeze-out ratio, but more importantly, of its sensitivity to the cosmic expansion. To 
cancel out the overall error of this model, only the difference X n (NNDE) — X n (ACDM) is 
used. 




-6 -5.9 -5.8 -5.7 -5.6 -5.5 -5.4 -5.3 



log 10 (m[eV]) 



Figure 5: type I NNDE constraint on g and the neutrino mass from the DE EOS and BBN 
measurements. The color-bar signifies the sigma confidence level. 

The 95% confidence level constraints of this model 2 , are therefore (see Figure 5), 

m v < 3-l(T 6 ey, 
#> 350. (4.5) 

2 The constraint on g is irrelevant, since it is a model related parameter. Using a different model, would 
have resulted in different constraints on its parameters, however, it would not have a significant effect on 
the neutrino mass constraint. 
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Figure 6: type II NNDE constraint on g and the lightest neutrino mass from the DE EOS and 
BBN measurements. The color-bar signifies the sigma confidence level. In this case, there is no 
constraint on the mass of the lightest neutrino. 



This is of the same order of magnitude as the value obtained in [17] It is worth noting that 
according to these constraints on the lightest neutrino mass and using the best estimate we 
have for the neutrino square mass difference (obtained from neutrino oscillations [10]), the 
muon neutrino mass is vn vyb ~ 8.7 • 10 _3 ey. Therefore, in accordance with the previously 
noted fact, that a particle does not contribute to DE, once the universe cools beyond 
a certain temperature T oc m, and knowing that the lightest neutrino dominates comic 
expansion at red-shifts z ~ 1 the muon-neutrino-induced DE, dominated the expansion of 
the universe, up to red-shifts of 

777/ 

* > ~ 3000 ~ z eq . (4.6) 

Here z eq is the red-shift of the radiation-matter equality. This is a remarkable result, 
because it tells us that the primordial DE would not have significantly affected structure 
growth, and yet leaves room for subtle affects, that will perhaps be measured in the future. 

4.3 Type II NNDE 

In this case, we require the electron-induced DE to be subdominant (compared to the 
radiation energy density) all the way up to the end of the freeze-out era T* ~ 200keV, so 

Cm\ « na rad T*\ (4.7) 
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with a ra d being the radiation density constant and k ~ 7.25 is the number of radiative 
degrees of freedom, at the corresponding temperature. Since, today we have Crri^ — /?de 5 
eq. 4.7 implies, for the neutrino mass, 




^ a rad 




4 • 10" 3 eK 



(4.8) 



For this mass range, we can no longer assume only the lightest neutrino contributes. There- 
fore, for this model, the calculations include all three neutrino generations. The relations 
between the different neutrino masses, were taken to be the nominal values from [10], 



Figure 6, presents the constraints of this model on the lightest neutrino mass and the 
parameter g. There are no constraints on the neutrino mass for this model. However, we 
see that the maximal allowed value of g is inversely proportional to m VT . In this case, 
all three neutrino species, are involved in generating todays DE. Therefore, there are no 
available particles to induce primordial DE in proximity to the radiation-matter equality, 
and thus, no observable effects on the cosmic structure formation. 

5. Summary and Discussion 

A natural model for neutrino dark energy (NNDE), free of exotic particles or couplings, was 
presented. Basing the theory only on naturalness and the assumption that todays DE is 
driven by the neutrino, several unique predictions were derived. A specific relation w a (wo) 
governing the linear evolution approximation of the EOS, follows from NNDE. Verifying the 
above relation, can strongly support this model. Although, quintessence thawing models 
can produce similar w^wo) dependance, distinguishing between the two may be possible 
via a more advanced analysis of observations, using second order expansion of the EOS. 

NNDE also requires the existence of primordial DE, driven by heavier particles, domi- 
nating the expansion of the early universe. To prevent the electron-induced primordial DE 
from affecting the BBN, two methods were proposed. The first, entitled type I NNDE, is to 
assure that the electron-induced DE decays prior to the BBN freeze-out phase. The second 
method, entitled type II NNDE, is to assure that the magnitude of the electron-induced 
DE, is small enough to have little affect of the BBN. Using the above methods, constrains 
on NNDE are placed. 

Type I NNDE, provides new constraints on the mass of the lightest neutrino, requiring 
m v < 3 • 10~ 6 eV. In addition, type I NNDE, predicts relic primordial DE, induced by the 
muon neutrino, in proximity to the radiation-matter equality phase. This would mildly 
affect structure formation and the CMB. Observing these effects could prove as the main 
finger-print of the NNDE. 



Aml e = 7.6 • 10~ 5 eV 2 , 
Ami,, = 2.4 • 10~ 3 eV 2 . 



(4.9) 
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